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A previous method for handling bound states in QCD is briefly revisited. Taking 
advantage of the Feynman-Schwinger representation for the iterated quark prop- 
agator in an external field, it is possible to give closed representations for certain 
appropriate (second order) two point and four point Green functions, H^^^(x — y) 
and (x\ ,X2,yi,y2), as path integrals on quark world lines. Then, starting from 
reasonable assumptions on the Wilson line correlators, a Bethe-Salpeter equation 
for iJ^*' and a Dyson-Schwinger equation for H^^^ can be obtained, which are 
consistent with the Goldstone theorem in the chiral limit. Such equations are 
too complicate to be solved directly. However, a reduced Salpeter equation can 
be derived which is tractable and has been applied to a calculation of the meson 
spectrum. The results are in general good agreement with the data, but with the 
important exceptions of the light pseudo scalars (that are related to the breaking 
of the chiral symmetry). In this scenario two important improvements can be in- 
troduced: a) the fixed coupling constant can be replaced by a running coupling 
constant as{Q'^) appropriately modified in the infrared region; b) the fixed mass 
in the reduced equation can be replaced for light quarks by an effective mass de- 
pending on the momentum of the particle, as suggested by the form of the DS 
equation. Then even the light pseudo scalar mesons can be made to agree with to 
their experimental value. 



1. Introduction 

In previous papers we have introduced wliat we liave called a second order 
Bcthe-Salpeter formalism, which works in terms of certain appropriate four 
point and two point Green functions, H^'^\xx, X2, 1/1,1/2) and H^'^\x ~ y)^ . 
Taking advantage of a Feynman-Schwinger representation for the "sec- 
ond order" quark propagator in an external field, it is possible to write 
i?^^^ and iJ*-^-* as path integrals on quark or antiquark world lines join- 
ing y to X, to xi_2. In such representations the gauge field appears 
only trough Wilson line correlators like ■^Trcoior(exp[i dz^^Af^{z)] or 
|Trcoior(exp[z dz'^A^lz)] exp[i JJ^^ dz^ A^{z)\) which reduce to ordinary 
Wilson loops W in the limit x ^ y or xi X2 and yi y2- In analogy 
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with a usual assumption on W such correlators are written as the sum of 
their perturbative expressions and appropriate area terms and it is possible 
to obtain a confining Bethe-Salpeter equation for the quantity H^^^ and a 
corresponding Dyson- Schwinger equation for H^^\ 

The above equations are too complicate to be solved directly. However, 
by a conventional three dimensional reduction, one can derive a tractable 
squared mass operator that can be applied to an evaluation of the spec- 
trum ^. With an appropriate choice of the parameters (quark masses, cou- 
pling constant, string tension) the results can be made in good agreement 
with the data in the entire framework of the light-light, light-heavy and 
heavy- heavy quark- antiquark sectors, with the important exception, how- 
ever, of the light pseudo scalar mesons ^. 

In this paper we want to discuss two important improvements to the 
above scenario: 

a) the fixed coupling constant as is replaced by a running coupling con- 
stant as{Q^), which amounts to take into account higher order contributions 
in as (and in particular closed quark loops); 

b) an effective mass for the light quarks is introduced which depends on 
the momentum of the particle, as suggested by the form of the DS equation. 

As well known, the usual perturbative expressions for as{Q^) have an 
unphysical singularity for = Aqqj-,, which would be disastrous for our 
purpose. However, various modification have been proposed for the in- 
frared region We have considered in particular the Shirkov-Solovtsov 
prescription, which rests only on general analyticity requirements. We find 
that with such prescription, and paramctrazing the effective squared light 
quark masses by a simple polynomial, even the light pseudo scalar mesons 
can be made in agreement with their experimental values. It should be 
stressed that in the fit only the quark masses are treated as free param- 
eters, while the constant A and a occurring in the BS kernel arc a priori 
fixed on the basis of high energy phenomenology and lattice simulations. 
Results in this line were already published in ref. [^]; here we present a 
more systematic study. 

The remaining part of the paper is organized as the following plan: in 
sections 2 and 3 we briefly revue the second order Bethe-Salpeter formalism, 
in section 4 we discuss the running coupling constant and efl[ective mass, in 
sections 5 and 6 we reports the results and make some conclusions. 
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2. Second order correlators 

After integrating out the fermionic fields, the appropriate ordinary (first 
order) four point function for the quark-antiquark bound state problem 
can be written 

G^^\xuX2;yi,y2) = ^ $^(0|TVla(a:i)?2a(a^2)?ifc(2/l)V'2h(2/2)|0) = 
ah 

= ~Trco\oT{Si{xi,yi;A)S2{y2,X2;A)), (1) 

Here a and h are color indexes, the subscripts 1 and 2 refer to the quark 
and the antiquark respectively, projection on the color singlet has been 
performed, S{x, y; A) denotes the quark propagator in an external field, 

{i'r^'D^-m)S{x,y;A) = 6\x-y), (2) 

and 

{/[A]) = J DAMf[A] e^^'^^^^flA] , (3) 

with 

Mp[A] = Detn2^i[l + gj^'A^iijPj^ - mj)-^]. 
The second order four point function is defined by 

H^^\xi,x2;yi,y2) = -^TT,^ioA^i{xuyi;A)A^{y2,x2;A)), (4) 

where Ai{x, y. A) is the second order propagator that satisfies the second 
order differential equation 

{D^D^ + m^ - ^gc7^^-'F^^)A^x,y;A) = -5\x - y) , (5) 

(cr''^ = |[7'',7'']) and is related to S{x,y;A) by 

S{x, y- A) = {lYD^ + m)A'^(a;, y; A) . (6) 

The quantity i?'^' is related to G^^^ by an integro-differential operator, of 
the type [(?7f 9i^ + mi)(i72 92,y + TO2) + - ■ •] , which we do not need to specify 
in detail. The important fact is that the two functions are completely 
equivalent for the determination of the bound states since their Fourier 
transforms H'^'^^ and G^^-' have related analyticity properties and the same 
poles. 
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The advantage in considering second order quantities is that it is possible 
to write for A"^ (x,y; A) a generaUzed Feynman-Schwinger representation, 
i. 6. to solve eq. (5) in terms of a quark path integral 



^ pOO nX f-S ^ 

A''{x,y;A) = -- / ds I £'zexp[-i / dT-{m^ 
2 Jo Jy Jo 2 

S'o^eMig f dTZ^'A^{z)], (7) 
Jo 

where the world-line zi^ = z'^{t) connecting y to x is written in the four- 
dimensional language in terms of an additional parameter r, 5q = T exp — 

\ drai^^jg^ and SSi^" = dz^'Sz" - dz^Szi^ (the functional derivative 
being defined through an arbitrary deformation, 2; ^ 2; + ^2;, of the world- 
line). 

Replacing eq. (7) in eq. (4) a similar representation can be obtained 
for the 4-point function 

H^^\xi,X2;yi,y2) = {T;f dsi ds2 Dz, Dz2 

^ JO Jo Jyi Jy2 

exp { - ^ / '^^i^™! + -^i) ~ ^ / '^^^(m^ + i^)} (8) 

-S'o'S'o'Ti,oior{Pexp{ig d^f A^(2i)}P exp {ig / dz^ Af,{z2)}) , 

<J Jyi Jx2 

The interesting aspect of the above equation is that the gauge field appears 
in it only through the expectation value of the product of the two Wilson 
lines. 

Similarly the second order two point function (uncolored full quark prop- 
agator) H^'^\x — y) = -^T^t: color {S{x,y : A)) can be written 

1 /■CO j-X ■ 

H^^\x-y) = - ds Dzexp{-- dT{m'^ + z^)} 

^ Jo Jy 2 Jq 

;^5'o'Trcoior(Pexp{ig^ dz^" A^{z)}) . (9) 
Notice that, for large separations, we can also write in operatorial form 
H^^\xi,X2;yi,y2) = ^ ^(0|T.^i„(a;i)^2a(a;2)?^i6(yi)</'26(t/2)|0) 



3 , 

ab 



and 

H('\x-y) = j2mM^)i^a{ym, 
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(j){x) being defined by 

V'(a;) = {ij"Di, + m)(j){x) . 

The above forms are important for a definition of the relativistic wave 
function. 

3. Bethe-Salpeter and Dyson- Schwinger equations 

In the hmit X2 ^ Xi, y2 ^ yi oi y ^ x the two Wilson hnes occurring in 
(8) or the single line occurring in (9) close in a single Wilson loop F 

W = (iTrcoiorPexp [ ^dz^'A^iz)]} . (10) 

As in the previous papers we assume that in a first approximation i In 
can be written as the sum of its perturbative expression and an area term 
ilnW = (ilnW^)pert + fS. Then, at the lowest order in the couphng con- 
stant, we can assume 

i\nW =^g'^ j) dz^" j) dz^'D^^iz- z') + (11) 

a^dz^^ dz''5{z' - z°')\. - z'l J^' dx{l - [X^ + (1 - A)^]^}^ . 

Notice that the surface term in (11) is written as the algebraic sum of 
successive equal time straight strips (zt denotes the transversal compo- 
nent of z). For a flat loop or for other special geometries (e. g. for two 
quarks uniformly rotating around their fixed center of mass) this coincides 
obviously with the plane or the minimum surface delimited by T. That is 
not generally the case, indeed the right hand side of (11) usually depends 
on the reference frame. Since, however, in contrast e. g. with S'min, such 
quantity maintains many of the analytic properties of the original i In W, 
we shall assume (11) to be valid for an arbitrary loop in the center of mass 
reference frame. Actually we shall assume this even for X2 ^ xi, y2 ^ yi or 
y ^ X, in analogy with what happens in the pure perturbative case. In this 
way single perturbative and confinement contributions are put on the same 
foot and we may refer to them as a gluon exchange and a string connection 
between the two quarks. 

Replacing (11) in (8) and (9) we obtain the following equations 

poo />oo />a:i />X2 

H^^\xi,X2;yi,y2) = i7;f dsi ds2 Dzi Dz2 

^ Jo Jo Jyi Jy2 
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and 



{ - ^ E I ' "^^^ + } ^o' ^o' E y/ '^^^ (12) 

y dT!jE{zj - Zj;Zj, z'j) -i J dn J dT2E{zi - Z2; ii, i2)| 

H^^\x-y) = 7; / ds / Dzexp{-'- / dT(m2 + i^)} 

^ Jo Jy ^ Jo 

5o^exp{i / / E{z-z';z,z')) , (13) 
JO Jo 

where we have set 

E{c,p,p') = £;pert(C;p,p') + -Bco„f(C;p,p') (14) 

with 

('i;pert = 4^|as-DM.(C)pV" . . 

\ i^conf = <5(Co)|Ck(Po)e(p^) /o rfA{pgp^,2 _ [^p/ + (1 _ A)popJr]'}' • ^ ^ 

From cqs. (12) and (13), by various manipulation and using an ap- 
propriate iterative procedure, a Bethe-Salpeter equation for the function 
if'^' {xi ,X2;yi, 2/2) and a Dyson-Schwinger equation for iJ^^^ (a; — y) can be 
derived in a kind of generalized ladder and rainbow approximation respec- 
tively in the form 

H^^'' {x, , X2 ; 2/1 , 2/2) = (xi - 2/1 ) i/f ^ (X2 - 2/2) - 

- i J d%d%d%d%H[^\xi - a) fff ^(X2 - 6) 

X lab (a , 6 ; »?i , ??2 ) 0-5* F^^'H'?! > ??2 ; 2/1 , 2/2) , (16) 

and 

F(2)(a; - 2/) = if,$''(a; - 2/) + » 1 d^d^r^dH' d'r,' H^^\x - 

x/a6(C,^';'?,V)«^"if^'n^ - »?')^'if'''(^' - 2/) , (17) 
where a, 6 = 0, /it', we have set cr° = 1, and and if2^' denote the 

second order full quark and the antiquark propagators respectively. 

In momentum representation, the corresponding homogeneous BS- 
equation can be written 

^(2) ( Ip + fc)a«$p(u)a*^f ( _ ip + fc) , (18) 
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the center of mass frame has to be understood, P = (mB,0) and ^p{k) 
denotes the appropriate second order wave function 

Similarly, in terms of the irreducible self-energy, defined by 

— k^ — 

the DS-equation can be written also 

m= j (0J-J4k-I;^,^)a^m^\l)a'> . (19) 

The kernels in (18) and (19) are the same in the two equations, consistently 
with the requirement of chiral symmetry limit ^, and are given by 

ioAQ;P,p') = 4 J d\e'QiE{C;p,p') = 167r^a,pV''^a/3(Q) + 

+4a J d^Ce-''^<\C\e{po)e{p'o) ^ dX{plp'o^ - [Xp'^PT + (1 - A)poPt]'}^ 
V;o(g;P,p') = 47r«^«s(<5;^0. - S^Q^)p'f,Do.0iQ) - 
-a [d'Ce-^^<e{po)^f^^=S^p'o 

io-MQlP,p') = -47^^^a,K('5;^Q. - S^Qp)D^(}iQ) + 



, f ,3/- -iQ-C ipPcr CcrPp , 



Wa(Q;p,p') = TT^Md^Q, - d^Q^^WpQ, - 5'^Qp)D^p{Q) , (20) 

where in the second and in the third equation Co = has to be understood. 

To find the qq spectrum, in principle one should solve first (19) and 
use the resulting propagator in (18). In practice this turns out to be a 
difficult task and one has to resort to the three dimensional equation which 
can be obtained from (18) by the so called instantaneous approximation. 
This consists in replacing Hj (k) in (18) with the free quark propagator 
-r^^ — J- and the kernels lab with its so called instantaneous approximation 

k -m- 

7X*(k,k') \ 
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The reduced equation takes the form of the eigenvalue equation for a 
squared mass operator = Mq + U, with Mq = ^mf + + 
and 

(k|[/ik') = 1 /^±^i^^st(k,k'),/5+5,^^^ (21) 

(for an expHcit expression see ref. [®'^]). The quadratic form of the above 
equation obviously derives from the second order formalism we have used. 

Alternatively, in more usual terms, one can look for the eigenvalue of the 
mass operator or center of mass Hamiltonian -ffcM = M = Mq + V, with 

V defined by MqV + VMq + V'^ = U. Neglecting the term V"^, the Hncar 
potential V can be obtained from U simply by the kinematic replacement 
I {wi+w2)(wi+w^) _^ — ^ Such expression is particularly useful for 

a comparison with models based on potential. In particular in the static 
limit V reduces to the Cornell potential 

V = -i^+ar, (22) 
6 r 

in the semirelativistic limit (up to ^ terms after an appropriate Foldy- 
Wouthuysen transformation) equals the potential discussed in ref. p]. If 

the spin dependent terms are neglected but full rclativistic kinematics is 
kept, it becomes identical to the potential of the (rclativistic) flux tube 
model ^. 

4. Running coupling constant and effective mass 

In ref. ["^J the spectrum was evaluated for both the operators and Hcu 
introduced in the preceding section, omitting the spin-orbit terms in the 
potential but including the hyperfine ones. As we told, with fixed coupling 
constant and quark masses, a general good fit of the data was obtainc;d over 
the entire calculable spectrum with, however, the relevant exception of the 
light pseudoscalar mesons. At the light of the idea that such mesons should 
be Goldston massless particles in the chiral limit, this is not surprising 
The results obtained in ref. [^] suggest, however, that the situation can be 
greatly improved using a running coupling constant and an effective light 
quark mass function of the momentum, as implied in (19) for a kernel of 
the type (20). 

At one loop the running coupling constant in QCD is usually written as 
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Q being the relevant energy scale, /?o = 11 — f^fi ^-i^d A^f the number of 
flavors with masses smaller than Q. Such expression becomes singular and 
completely inadequate as approaches A^. 

The most naive modification of eq. (23) consists in saturating a^iQ'^) 
to a certain maximum value ois as decrease (fig. 1) and in treating this 
value as a phenomenological parameter (truncation prescription). 

Alternatively Shirkov and Solovtsov ^ replace (23) with 



Po 



1 



A2 



(24) 



This remains regular for 
a,{0) = 47r//3o, for ^ 0. 



Jn(Q2/A^) A2-Q2^ 
= A^ and has a finite A independent limit, 
Eq. (24) is obtained assuming a dispersion 
relation with a cut from — oo < < and applying (23) to the evaluation 
of the spectral function alone. 

In the quark-antiquark bound state problem the variable can be 
identified with the squared momentum transfer = (k — k')^. Then, 
typically, (Q^) ranges between (0.1 GeV)^ and (1 GeV)^ for different 
quark masses and internal excitations and values of smaller than A^ 
can be important. The specific infrared behavior is therefore expected to 
affect the spectrum and other properties of mesons. 




Q (GeV) 



Figure 1. Running coupling constant as(Q^) on logarithmic scale. Perturbative expres- 
sion (dotted line), truncation prescription (dashed line), Shirkov-Solovtsov prescription 
(full line), Nesterenko prescription (dot-dashed line). 



Coming to the quark masses, notice that, even neglecting the spin de- 
pendent part in H^'^\k) (consistently with what we do for (k|J7|k')), we 
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could more sensibly write H^'^^k) = i/{k^ — — ro(fc)) in the reduction 
procedure, rather than simply identify such quantity with its free expres- 
sion (ro(fc) = 0). However, since the Ia{Q,p,p') are not formally covariant, 
but given in terms of CM. variables, ro(fc) must depend separately on 
and k^. Then the pole of i?2(fc), defined by 

A;g-k2-m2-r(A;g,k2) =0, (25) 

could be written 

fco=meff(|k|) + k^ (26) 

where mgg(|k|) would be a |k| dependent expression that is expected to 
approach the current for |k| — > and to increase toward a kind of 
constituent m'^ as |k| increases Then, eventually, we obtain the same 
operator as given above, but with mi and m2 replaced by expressions 
like mlM- 

As matter of fact we have not tried in this paper to evaluate r(p) from 
first principles by solving eq. (19) but have simply parametrized mgjj(|k|) 
as a polynomial in jk| for the light quarks and used a fixed mass as in the 
preceding papers for strange and heavy quarks. 

5. Numerical procedure and results 

In the calculation of the eigenvalues of the operator introduced in sect. 
4 we have used the same approximations and the same numerical procedure 
as in the previous papers. Wc have retained only the hypcrfinc part of the 
complicate spin dependent terms (spin orbit and tensorial ones) occurring 
in the potential U. We have solved the eigenvalue equation for Hqm with 
V equal its static limit (22) by the Rayleigh-Ritz method, using the three 
dimensional Harmonic oscillator basis and diagonalizing a 30 x 30 matrix. 
Then we have evaluated {4>n\M'^\(j>n) for the eigenfunction obtained in 
the first step. 

In figs. 2 and 3 wc have represented graphically the results already 
obtained in ref. p] with a truncated running coupling constant (crosses) and 
the new ones performed with the Shirkov-Solovtsov prescription (circlets), 
contrasted with the data (small lines) [^]. In tables I- VII we give the 
corresponding numerical values. 



''Take into account that from the virial theorem for a hnear potential in the extreme 
relativistic approximation we have (|k|) = a{r) and a large |k| corresponds to a peripheral 
interaction or, what is the same, to a small |Q|. 
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(GeV) 




Figure 2. Light-light quarkonium spectrum and Rcggc trajectories, rtj = 4, A = 
0.2 GeV. Crosses and dashed line: truncation prescription, <5s = 0.35, cr = 0.2 GcV^, 
ms = 0.2 GeV, ruu = rud = 0.01 GeV. Circlets and full line: Shirkov-Solovtsov as{Q^), 
c = 0.18 GeV^, ms = 0.381 GeV, = m| = 0.121 k - 0.025 k"^ + 0.25 k'^. 



In both cases we have taken TVf = 4 and A = 200 MeV. With the 
truncation prescription, the Hght quark masses were fixed on typical CTirrent 
values, mu = Wd = 10 MeV, = 200 MeV, the maximum value for 
Q!s(<3^)j the string tension and the heavy quark masses were treated as 
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(GeV) 



^So ^Sl 



n 

> 6 

> 5 
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'1 -fj 

n 
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bb 



(GeV) 



^So ^Si P 


^So^Si P 
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n n 


2. 


- o 2 " ° 

1 


qb 
q = u,d 


sb 



iSo ^Si iPi ^Pj 



n 



n n 

J J 



° 2 - 
° 1 : 



, .2 

1 1 





cc 





p 


^So^Si P 


n n 




n n 






> o2" ° ^ 








q = 


u, d 


sc 



Figure 3. Hcavy-licavy and licavy-light quarkonium spectrum. Uf = A, A = 0.2 GeV. 
Crosses: truncation prescription, <5s = 0.35, cr = 0.2 GcV^, mj, = 4.763 GcV, rric = 
1.394 GcV, nis = 0.2 GeV, m„ = = 0.01 GeV. Circlets: Shirkov-Solovtsov a^Q^), 
(T = 0.18 GeV^, mi = 4.898, rUc = 1.545, = 0.381 GeV, ml = ml = 0.121 fc - 
0.025 fe2 + o.25 k*. 



adjustable parameters and chosen as as = 0.35, cr = 0.2 GeV^ (in order to 
reproduce the correct J/\l/ — rjc separation and Regge trajectory slope) and 
rric = 1.394 GeV, rrib = 4.763 GeV (in order to obtain exact masses for J/^ 
and T). On the contrary, with the Shirkov-Solovtsov prescription, wc have 
fixed the string tension on the commonly accepted value cr = 0.18 GeV^ 
and used all masses as adjustable parameters. For the light quarks we have 
taken (in GeV^) 

ml = ml = 0.121 k - 0.025 + 0.25 k^ , (27) 

(fc being the quark momentum in the meson C. M. frame), for the 
strange and the heavy quarks = 0.381 GeV, rric = 1-545 GeV and 
nib = 4.898 GcV. 

As it is apparent the use of a running coupling constant with infrared 
truncation does not improve essentially the situation in comparison with 
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Table I. qq {q = u,d), Uf — 4, A ^ 0.2 GcV. (a) truncation prescription 
tts = 0.35, cr = 0.2 GeV^ m„ = = 0.01 GeV. 





(b) Shirkov-Solovtsov Q!s((3^), a = 


0.18 GeV , 






2 2 


= 0.121 k — 0.025 fc + 


r\ 7,4 




States 




experimental values 


(a) 


(b) 






(MeV) 


(McV) 


(MeV) 


l^So 




134.9764 ± 0.0006 


479 

1 


135 




139.56995 ± 0.00035 _ 


l^Si 


p(770) 


768.5 ± 0.6 


846 


739 


lASS 




630 


367 


604 


2^80 


7r(1300) 


1300 ± 100 


1326 


1402 


2 3Si 


p(1450) 


1465 ± 25 


1461 


1509 


2ASS 




165 


135 


107 


S^So 


7r(1800) 


1795 ± 10 


1815 


1994 


S^Si 


p(2150) 


2149 ± 17 


1916 


2064 


3ASS 




354 


101 


70 


l^Pi 


61(1235) 


1231 ± 10 


1333 


1310 


1^P2 


02(1320) 


1318.1 ±0.7^ 








ai(1260) 


1230 ±40 > 1303 


1333 


1324 


l3Po 


0,0 (1450) 


1450 ± 40 J 








7r2(1670) 


1670 ± 20 


1701 


1739 




P3(1690) 


1691 ±5 ^ 






13D2 






1701 


1743 


l^Di 


p(1700) 


1700 ± 20 J 






l^Fa 






1990 


2043 


13F4 


a4(2040) 










X(2000) 


2037 ± 26 1 


1990 


2044 












IIG4 






2238 


2318 


l^Gs 


P5(2350) 








13G4 




2330 ± 35 1 


2238 


2321 


l^Gs 


/93(2250) 














2460 


2570 




06(2450) 












2450 ± 130 1 


2460 


2570 


13H4 
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Table II. ss, Uf = 4, A — 0.2 GcV. (a) truncation prescription ag = 0.35, 
a = 0.2 GeV^ = 0.2 GeV. (b) Shirkov-Solovtsov a^iQ^), 
C7 = 0.18 GeV^ rris = 0.381 GeV. 



States 




experimental values 


(a) 


(b) 






(MeV 




(MeV) 


(MeV) 


1 bo 


Vs 


r r;(547) ] 
I V(958) J 


•674 


/ yu 






(^(1020) 


1019.413 ± 0.008 


1049 


1010 


lASS 




335.3 ± 0.1 


259 


313 










1503 


1469 


2% 


(j!)(1680) 


1680 ± 20 


1630 


1600 


2ASS 








127 


131 


3IS0 








1977 


1934 


3% 








2068 


2023 


3ASS 








91 


89 


l^Pi 


/ii(1380) 


1380 ± 20 


1491 


1484 


13P2 


/2(1525) 


1525 ± 5 ' 








l3Pi 


/i(1510) 


1512 ±4 


> 1518 


1492 


1496 


l^Po 


/o(1500) 


1503 ± 11 ^ 








1^D2 


r?2(1870) 


1881 ± 32 ± 40 


1830 


1826 


l^Dg 


03(1850) 


1854± 7^ 






l^Da 








1830 


1829 


l^Di 






1 






l^Fs 








2103 


2089 


13F4 


/j(2220) 


2225 ± 6 






l^Fg 






1 


2103 


2091 


l^Fa 


/2(2150) 










IIG4 








2339 


2313 


l^Gs 












13G4 




1 




2339 


2316 


l^Gg 












l^Hs 








2553 


2516 


















) 




2553 


2516 


13H4 













the fixed coupling constant case. The same happen with other prescrip- 
tion like the Dokshitzer et al. prescription ^. On the contrary with the 
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Table III. qs {q = u,d), Uf =4, A = 0.2 GcV. (a) truncation prescription 
as = 0.35, cr = 0.2 GeV^ = 0.2 GeV, to„ = = 0.01 GeV. 
(b) Shirkov-Solovtsov a^iQ'^), a = 0.18 GeV^ = 0.381 GeV, 





2 2 


= 0.121 k — 0.025 k + 


0.25 k . 




States 




experimental values 


(a) 


(b) 






(MeV) 


(MeV) 


(MeV) 


1 ^Sn 




497.672 ±0.0311 
493.677 ±0.016 J 


645 


470 


l^Si 


if* (892) 


891.59 ± 0.24 


943 


868 


lASS 




393.9 ± 0.3 


298 


398 


2^80 


K{U60) 




1411 


1426 


23Si 


K*{U10) 


1412 ± 12 


1541 


1545 


2ASS 






130 


119 


3IS0 






1892 


1956 


S^Si 






1989 


2035 


3ASS 






97 


80 


l^Pi 


iri(1270) 


1273 ± 7 


1412 


1397 


13P2 


/r|(1430) 


1425.4 ± 1.3 1 






l^Pi 


^^1(1400) 


1402 ±7 > 1418 


1412 


1410 


l^Po 


i4:o*(1430) 


1429 ±5 J 






l^Da 


7^2(1580) 




1766 


1775 


l^Ds 




1770 ± 10 ^ 






13D2 


if 2 (1770) 


1773 ± 8 i 1765 


1766 


1779 


l^Di 


if* (1680) 


1714 ±20 J 






l^Fs 






2046 


2067 


13F4 


if|(2045) 


2045 ±9 ^ 






l^Fg 


is:3(2320) 


2324 ±24 > 2130 


2046 


2069 


l^Fa 


i4:|(1980) 


1975 ±22 J 






IIG4 






2289 


2317 


l^Gs 


if5*(2380) 


2382 ±14 ±19^ 






1^64 


if 4 (2500) 


2490 ± 20 > 


2289 


2320 


l^Ga 










l^Hs 






2506 


2545 


l^He 










l^Hs 




1 


2506 


2545 


13H4 
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Table IV. 66, Uf — A, A = 0.2 GeV. (a) truncation prescription as = 0.35, 
cr = 0.2 GeV^ rub = 4.763 GeV. (b) Shirkov-Solovtsov a^iQ'^), 
a = 0.18 GeV^ = 4.898 GeV. 



St at PR 




PXTiPriTTifnt^il VriliiPS 


Ca) 


fbl 






(MeV) 


(MeV) 


(MeV) 


l^So 






9374 


9374 


l^Si 


T(15) 


9460.30 ± 0.26 


9460 


9460 


2^So 






9975 


9988 


2 3Si 


T(25) 


10023.26 ± 0.31 


10010 


10023 


3 ^So 






10322 


10342 


3 3Si 


T(3S') 


10355.2 ± 0.5 


10348 


10368 


4IS0 






10598 


10618 


4% 


T(45) 


10580.0 ± 3.5 


10620 


10639 


5IS0 






10837 


10854 




T(10860) 


10865 ± 8 


10857 


10872 


6^80 






11060 


11070 


6 3Si 


T(11020) 


11019 ± 8 


11079 


11089 


l^Pi 






9908 


9918 


13P2 


Xb2{lP) 


9912.6 ±0.5^ 






l3Pi 


XbiilP) 


9892.7 ±0.6 \ 9900 


9908 


9920 


l^Po 


XboilP) 


9859.9 ± 1.0 J 






2iPi 






10260 


10279 


2 3P2 


Xb2{2P) 


10268.5 ±0.4^ 






2 3Pi 


xbiim 


10255.2 ±0.5 \ 10260 


10260 


10280 


2 3Po 


X6o(2P) 


10232.1 ±0.6 J 







Shirkov-Solovtsov coupling constant, together with the running u, d mass 
of eq. (27), the main difficulties seem to be solved. Even the masses of 
TT, r]s and K turn correct without destroying the agreement in the other 

part of the spectrum. An exception could be represented by the masses of 
the uh channel that for some reason turn out to be high. The use of an 
effective running mass even for does not seem to produce any substan- 
tial modification. With reference to the mass of r?s reported as a data in 
fig. 2 concerning the ss sector, it should be mentioned that the quantity is 
derived from the experimental values of rj and 77' masses with conventional 
assumptions on the origin of this particles as a mixing of 4= [uu + dd) and 
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Tabic Y. cc, Uf = 4, A — 0.2 GcV. (a) truncation prescription = 0.35, 
cr = 0.2 GeV^ mc = 1.394 GeV. (b) Shirkov-Solovtsov a^iQ^), 
a = 0.18 GeV^ rUc = 1.545 GeV. 



States 




experimental values 


(a) 


(b) 






(MeV) 


(MeV) 


(MeV) 


l^So 


^c(15) 


2979.7 ± 1.5 


2982 


2977 


1% 


J/^(15) 


3096.87 ± 0.04 


3097 


3097 


lASS 




117 


115 


119 


2^80 


r?e(25) 


3594 ± 5 


3575 


3606 


2 3Si 


^(25) 


3685.96 ± 0.09 


3642 


3670 


2ASS 




92 


67 


64 


3 ^So 






3974 


4005 


33Si 


^/'(4040) 


4040 ± 10 


4025 


4054 


4% 






4298 


4323 


4% 


^(4415) 


4415 ± 6 


4341 


4364 


l^Pi 






3529 


3556 


l^Ps 


Xc2(lP) 


3556.18 ±0.13^ 






l3Pi 


Xciim 


3510.51 ±0.12 I 3525 


3530 


3561 


l^Po 


Xco(lP) 


3415.1 ±0.8 J 






2iPi 






3925 


3954 


2 3p 






3927 


3958 


l^Da 






3813 


3853 












l^Da 


V'(3836) 


3836 ± 13 I 


3813 


3854 


l^Di 


V'(3770) 


3769.9 ± 2.5 J 






2ID2 






4149 


4183 


2^03 










2^D2 




4159 ±20} 


4149 


4184 


2 3Di 


^(4160) 









ss states. 

6. Conclusions 

In conclusion, starting from a reasonable ansatz on the Wilson loop it is 
possible to develop a second order Bethe-Salpeter approach to the qq bound 
state problem. From this, by a standard reduction technique, an eigenvalue 
equation for a mass squared operator M"^ = {w\ + w^f + U is obtained, 
that can be effectively applied to the calculation of the spectrum. 
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Table VI. Light-heavy quarkonium systems, n/ = 4, A = 0.2 GeV. 
(a) truncation prescription as = 0.35, a = 0.2 GeV^, rrib = 4.763 GeV, 
rUc = 1.394 GeV, m„ = = 0.01 GeV. (b) Shirkov-Solovtsov as(Q^), 
cr = 0.18 GeV^ rrib = 4.898 GeV, = 1.545 GeV, 
ml = ml = 0.121 k - 0.025 P + 0.25 k'^. 



States 




oxpcrimcntal vahics (a) 


(b) 






(McY) 




(McV) 


qc {q = u, d) 










1 la 
1 So 




1869.3 ±0.5' 
1864.5 ±0.5^ 


> 1875 


1862 


1 3c 
1 ''bi 


f L»*(2010)± 
D*(2007)° 


2010.0 ±0.5' 
2006.7 ±0.5^ 


>■ Z\JZ\J 


on 1 

zOlz 


lASS 




141 ± 1 


145 


150 


2^80 


£)' 


2580 


2525 


2615 


23Si 


D*' 


2637 ± 8 


2606 


2683 


2ASS 




57 


81 


67 


l^Pi 






2474 


2504 


13P2 


f £)*(2460)± 


2459 ± 4 ^ 






D^(2460)° 


2458.9 ± 2.0 






l^Pi 


'Di(2420)± 


2427 ± 5 


2475 


2511 


_ 1)1(2420)° 


2422.2 ± 1.8 






l^^Po 




> 






(76 (<7 = u, d) 










l^So 




5278.9 ± 1.8' 
5279.2 ± 1.8^ 


> 5273 


5315 


l^Si 


B* 


5324.8 ± 1.8 


5339 


5376 


lASS 




46 ± 3 


66 


61 


2^80 






5893 


6034 


23Si 


B*' 


5906 ± 14 


5933 


6064 


l^Pi 






5791 


5838 


13P 




5825 ± 14 


5792 


5841 



If only the one gluon exchange contribution is included in the pertur- 
bative part of the Bethe-Salpctcr kernel, but a truncated running constant 
Q!s((5^) is used, a reasonable good reproduction of the data is already ob- 
tained (as far as calculable) with the relevant exception of the light pseudo 
scalar mesons. However if the 8hirkov-8olovtsov prescription for as(Q^) 
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Table VII. Light-heavy quarkonium systems, n/ = 4, A = 0.2 GeV. 
(a) truncation prescription ag = 0.35, a = 0.2 GeV^, nib = 4.763 GeV, 

rUc = 1.394 GeV, = 0.2 GeV. (b) Shirkov-Solovtsov as(Q^), 
cr = 0.18 GeV^ rub = 4.898 GeV, = 1.545 GeV, nis = 0.381 GeV. 
States experimental values (a) (b) 



(MeV) (MeV) (MeV) 



sc 










I'So 




1968.5 ± 0.6 


1982 


1953 


l^Si 




2112.4 ± 0.7 


2120 


2102 


lASS 




144 ± 1 


138 


149 


2^80 






2617 


2610 


23Si 






2698 


2691 


l^Pi 






2547 


2555 


13P2 


Dsj(2573)± 


2573.5 ±1.7 ' 






1-^Pl 


£)si(2536)± 


2535.35 ±0.34 


> 2548 


2562 


l^Po 










sb 










l^So 




5369.3 ± 2.0 


5364 


5348 


l^Si 


b: 


5416.3 ± 3.3 


5429 


5416 


lASS 




47 ± 4 


65 


68 


2^80 






5985 


5969 


2^Si 






6024 


6009 


l^Pi 






5858 


5857 


13P 




5853 ± 15 


5859 


5860 



is adopted and a phenomenological running mass is used for the u and d 
quarks, even the light pscudo scalar mesons occur in the correct range. 

An open problem is a resolution of the complicated Dyson-Schwinger 
equation in some reasonable approximation (19) and an actual derivation 
of an equation of the type (27) from first principle. 

References 

1. N. Brambilla, E. Montaldi, G.M. Prosperi, Phys. Rev. D 54 (1996) 3506; 
G.M. Prosperi, in Problems of Quantum Theory of Fields, Pag. 381, B.M. 
Barbashov, G.V. Efimov, A.V. Efremov Eds. JINR Dubna 1999, hep- 
ph/9906237. 

2. For a different approach: Yu.S. Kalashnikova, A.V. Nofediev, Yu.A. Simonov, 
Phys. Rev. D 64 (2001) 014037; Yu.S. Kalashnikova, A.V. Nefediev, Phys. 
Lett. B 530 (2002) 117, and references therein. 



February 1, 2008 22:18 WSPC/Trim Size: 9in x 6in for Proceedings 



prcgiapl 



20 

3. M. Baldicchi, G.M. Prosper!, Phys. Rev. D 62 (2000) 114024; Fizika B 8 
(1999) 2, 251; M. Baldicchi, in "QCD: Perturbative or Nonperturbative ?" 
Pag. 325, L. S. Ferrcira, P. Nogueira, J. I. Silvan-Marcos Eds. World Scientific 

(2001) , hep-ph/9911268. 

4. A.L Sanda Phys. Rev. Lett. 42 (1979) 1658; Yu.L. Dokshitzer, A. Lucenti, 
G. Marchesini, G.P. Salam, JHEP 9805 (1998) 003; Yu.L. Dokshitzer, in 
29*'' International conference on High-Energy Physics (ICHE 98), Vancouver, 
Canada, A. Astbury, D. Axon, J Robinson Eds. World Scientific (1999), hep- 
ph/9812252; Yu.L. Dokshitzer, V.A. Khoze, S.I. Troyan, Phys. Rev. D 53 
(1996) 89. 

5. D.V. Shirkov, I.L. Solovtsov, Phys. Rev. Lett. 79 (1997) 1209; Theor. Math. 
Phys. 120 (1999) 1220; K.A. Milton, I.L. Solovtsov, O.P. Solovtsova, Phys. 
Rev. D 64 (2001) 016005; N.G.Stcfanis, W. Schrocrs, H.-Ch. Kim, Phys. 
Lett. B 449 (1999) 299, see also Eur. Phys. J. C 18 (2000) 137; A. V. 
Nesterenko, Phys. Rev. D 62 (2000) 094028; D 64 (2001) 116009; A. V. 
Ncstcrcnko, Mod. Phys. Lett. A 15 (2000) 2401; A. V. Nesterenko, I.L. 
Solovtsov, ibid. A 16 (2001) 2517; S. J. Brodsky, S. Mcnkc, C. Merino, 
Phys. Rev. D 67 (2003), 055008. 

6. M. Baldicchi and G. M. Prosperi Phys. Rev. D 66 (2002) 074008. 

7. M.B. Hecht, CD. Roberts, S.M. Schmidt, Phys. Rev. C 63 (2001) 025213, 
and references therein. 

8. A. Barchielfi, E. Montaldi, G.M. Prosperi, JVucI Phys. B 296 (1988) 625; 
Erratum-ibid. B 303 (1988) 752; A. BarchieUi, N. Brambilla, G.M. Prosperi, 
IlNuovo Cimento 103 A (1990) 59; N. Brambilla, P. ConsoH, G.M. Prosperi, 
Phys. Rev. D 50 (1994) 5878. 

9. The Review of Particle Physics, K. Hagiwara et al, Phys. Rev. D 66 010001 

(2002) . 



